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Abstract 

We use a diagrammatic approach to study color-neutral heavy particle production in 
nucleus-nucleus collisions in a quasi-classical approximation without small- £ evolution. In 
order to treat the two nuclei symmetrically, we use the Coulomb gauge which gives the ap- 
propriate light cone gauge for each nucleus. The resulting cross section is factorized into a 
product of two Weizsacker- Williams gluon distributions of the two nuclei when the trans- 
verse momentum of the produced scalar particle is around the saturation momentum. We 
confirm our results in covariant gauge where the transverse momentum broadening of hard 
gluons can be described as a diffusion process. The transverse momentum factorization 
manifests itself in light cone gauge but not so clearly in covariant gauge. 



1 Introduction 



Transverse momentum dependent (TMD) parton distributions are necessary to study less inclu- 
sive processes in high energy heavy ion collisions [lj. TMD parton distributions contain more 
information about the parton small- x degrees of freedom than the integrated parton distribu- 
tions and require a new factorization formalism which is different from the traditional collinear 
factorization [2]. TMD's have attracted much interests in recent years [6j|26]. There are two 
different unintegrated gluon distributions usually encountered in the literature [HIE]. The con- 
ventional unintegrated gluon distribution is related to the color dipole scattering amplitude, 
and most known processes are related to it [5]. Another unintegrated gluon distribution is 
constructed by solving classical Yang-Mills equations of motion in the McLerran-Venugopalan 
model, a simple model used to study the small- x gluon distribution in the saturation region for 
large nuclei [Hill]- In a light cone gauge, the valence quarks of a nucleus are treated as static 
color sources, and the gluon fields of the sources add coherently then become the non-Abelian 
Weizsacker- Williams (WW) A^ w field of the nucleus. The WW gluon distribution is quite dif- 
ferent from the conventional unintegrated gluon and a few physical processes which can directly 
probe this gluon distribution are known [5]. 

In this paper, we study nucleus- nucleus collisions analytically. We always assume the 
strong coupling constant a < 1, but take a 2 ^ 1 / 3 , with A the atomic number of a nucleus, 
to be a large parameter. We further assume that there is no QCD evolution in the gluon 
distributions of nucleons in the nucleus. Instead of considering difficult processes, such as gluon 
production [7H9], whose final state interactions are very complicated, we will limit ourselves 
to a much simpler case where the produced particle is color-neutral. We always assume that 
the mass of the scalar particle, M, is much larger than the saturation momentum Q s . The 
effective interaction we will use is £ e // = —\g<j> ( t > F^ u F a ^ lJ , where is a scalar field. This 
effective interaction is used to studied Higgs production in pp and pA [T9 ||2TH23] . however, the 
phenomenological aspect is not our major interest in this paper. The advantage of using such 
an effective interaction in a study of nuclei collisions is that we do not have to worry about 
the final state interactions of the produced particle, all the effects are from the initial state of 
the nuclei. Therefore, we can view the process as follows: the gluon fields in the two nucleons 
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are developed independently and gluon distributions are formed within each nucleus before 
the collision. During the collision the gluon fields are released by the nuclei and the scalar 
particle is produced. After the collision the nuclei have to rebuild their gluon fields which 
would induce further gluon radiation but it is irrelevant to the current problem because the 
particle is color neutral and knows nothing about the rebuilding processes. This simple physical 
picture is exactly what we have in the Color Glass Condensate framework, and the transverse 
momentum distribution that manifests itself in this process is the WW gluon distribution, which 
comes from neglecting the rebuilding processes of the gluon field of the nucleus. In a light cone 
(LC) gauge , the cross section, Eq. is factorized into a production of two unintegrated 
gluon distributions Eq. f lTUj) and involves a new gluon distribution Eq. ffT9l) . the so-called linear 
polarized gluon distribution [20]. In different kinematic regimes the cross section takes different 
forms. When Z 2 ~ Q 2 , where [ is the transverse momentum of the scalar particle, the cross 
section becomes a product of two WW gluon distributions of the two nuclei in coordinate 
space Eq. (I2T]) . Therefore, we also manage to find a process that can probe the WW gluon 
distribution in addition to the dijet production discussed in [5]. If Z 2 ^> Q 2 , the cross section is 
factorizable as given by Eq. (|27j) . 

The major question in the study of nuclear collisions is deciding what gauge to use. Gluon 
distributions in the two fast moving nuclei cannot be described naturally in one single LC gauge. 
In order to treat the two nuclei symmetrically we use Coulomb gauge which is the essential 
ingredient of our calculation. The advantage of using Coulomb gauge is that the Coulomb 
gauge propagator connected to a object with a big plus momentum component is equivalent to 
the A + = gauge propagator. Similarly for a left moving object Coulomb gauge is equivalent 
to A_ = gauge [TH[T5]. It is not surprising that one can do the calculation in this mixture 
of gauges. After all the A^ w field has only transverse components and is obtained from the 
classical equation of motion in A_ = A + = gauge ^U\. However, there are some technical 
subtleties behind this naive connection which will be discussed in detail in Section [2j 

Instead of trying to solve the classical field equations of the two colliding nuclei, in Section 
[3]we will use a diagrammatic approach to study color-neutral particle production in the nucleus- 
nucleus collisions. In order to simplify the diagrams the Slavnov- Taylor- Ward (STW) identities 
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[T8] will be heavily used. The diagrammatic STW identities were used in various contexts, 
for example, in obtaining a non-linear gluon evolution equation [13], in studying the quantum 
structure of WW fields [11], and in a neutral current DIS off a large nucleus [12]. We will start 
with some low order examples, then illustrate how to generalize from only one nucleon in each 
nucleus to an arbitrary number of nucleons in each nucleus and how the contribution of all soft 
gluons are resummed by A^ w . The cross section is manifestly f actor izable in the light cone 
calculation. In Section H] we also perform a calculation in covariant gauge as a confirmation of 
the result we have in the LC calculation. The physical picture is slightly different in covariant 
gauge. A hard gluon released by a nucleus will be multiply scattered, both elastically and 
inelastically, by nucleons of the other nucleus as the two nuclei pass through each other. The 
transverse momentum distribution of the hard gluon will be gradually broadened by the multiple 
scattering and can be found by solving a diffusion equation [T21I21]- If we sum up all possible 
nucleons from which the hard gluon could come from in a nucleus the WW gluon distribution 
shows up right before the collision, however the cross section does not manifestly factorizable 
at first glance. This is additional evidence that the WW gluon distribution is the right type of 
gluon distribution for this process. Moreover, comparing the calculation from these two different 
approaches, we believe that transverse momentum factorization can be easily achieved in light 
cone gauge, or more precisely in Coulomb gauge. Therefore light cone gauge with appropriate 
regularization in the propagator and the Ward identities are the two essential ingredients in 
obtaining a transverse momentum factorized formula. To our knowledge this is the first process 
where an exact analytic formula has been formed for a physical process, involving momenta on 
the order of Q s , in nucleus-nucleus collisions in the quasi-classical approximation. We might 
hope that this method can be further applied to study more complicated processes, such as 
gluon production. 



2 Coulomb and Light Cone Gauge 

In the LC gauge A ■ n = 0, n 2 = 0, the Yang-Mills field propagator contains a singular 
denominator l/(k ■ n) and must be regularized. The propagator that will be used in the 
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following calculation takes the form 



—I 



n a kp npk, 



(1) 



k 2 + ie 



n ■ k + ie n ■ k — ie 



where the momentum k flows from a to ft. Note that, due to the opposite direction of the 
momentum flow, the second and the third terms in the propagator always have the opposite ie 
pole in the complex k ■ n-plane. This way of regularizing the propagator was previously derived 
in canonical quantization [28] and in path integral formalism [27] and more recently in a gauge 
transformation analysis [2S] • Moreover, it has been realized that the ie prescriptions in the LC 
propagator are closely related to initial and final state interactions [T2 | [T3 t [T7 ] [26] . That is the 
ie poles in the light cone propagator dictate how gluons propagate in the x_- or x + -directions 
when they are emitted from physical objects, like nucleons. For example, in = n ■ A = A + 
gauge, the propagator (OQ) tells us that the gluon fields coming from a right-moving nucleon 
only propagate in the positive ^--direction, or in terms of gauge rotation, gluons coming from a 
source with coordinate x°_ only gauge rotates other sources that satisfy x_ > x°_. With gluons 
only propagating in one certain direction the ie prescriptions can be used to avoid initial or 
final state effects. In a neutral current DIS off a large nucleus the regularization can be chosen 
such that it can eliminate final state interactions |12] . Implementing different prescriptions, 
one can also shift a final state interaction to an initial state interaction [26] . 

However, Eq. ([1]) cannot be used immediately in the calculation of nucleus-nucleus colli- 
sions, because one single LC gauge choice would make the two nuclei asymmetric. For example, 
if we choose A + = gauge, then a right moving nucleus having a big plus momentum compo- 
nent naturally fits in with such a gauge choice. However, for a left moving nucleus it is awkward 
to develop gluon fields in such a gauge. The same problem arises if we want to use A_ = 
gauge. Therefore one single LC gauge, either A + = or A_ = 0, will make the two fast moving 
nuclei asymmetrical. Naturally one wants to use A + = gauge for the right-moving nucleus 
while A- = gauge for the left-moving nucleus, so that gluon fields can be developed naturally 
in each nucleus. Indeed, the gauge choice satisfying this need exists, although it is not obvious 
at first sight. It is the Coulomb gauge V ■ A = [14], in which the propagator reads 




(2) 
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where N ■ v = Vq for any vector v. When this propagator is connected to a fast right moving 
object, so that 3> fc 2 , fc 2 , the dominant part in Eq. (F2J) is D_p = ikp /(k 2 k + ) which coincides 
with the one from the Eq. (JT]) in A + = gauge except for the ie pole. Similarly for a left moving 
system Coulomb gauge is equivalent to v4_ = gauge. Therefore, in Coulomb gauge the gluon 
fields emitted from a right moving system are in A + = gauge, while for a left moving system 
in A_ = gauge. 

However, the Coulomb gauge does not indicate an explicit choice of ie's in the propagator. 
So if we want to choose Coulomb gauge as our overall gauge choice we have to show that LC 
calculation is independent of the ie prescriptions in the propagator. Since we will use two 
different LC propagators for the two colliding nuclei, there are totally four different choices of 
ie's but only three of them are different. Different choices of ze's correspond to different ways 
that gluons propagate, i.e. different way of calculating diagrams. In principle we can compare 
diagrams coming from different ie choices order by order. Up to tree level diagrams, we manage 
to show that three different choices of ze's give an identical result. For the problem at hand we do 
not need to go beyond tree level calculation, since it is sufficient to compare tree level diagrams 
in the quasi-classical approach. Part of our goal here is quite similar to what has already 
been done by the authors in Ref. [26], where they introduced a gauge invariant transverse 
momentum-dependent parton distribution and showed that the definition is independent of the 
ie prescriptions on the LC propagator and is the same as the calculation done in covariant gauge. 
Here, however, we are not only able to show that WW gluon distribution is independent of the 
ie prescriptions, but also that the cross section for scalar particle production is independent of 
the prescriptions and factorized. 



3 Light cone calculation of scalar particle production 

In this section we evaluate scalar particle production in a nucleus-nucleus collision in light cone 
gauge. In light cone gauge noncausal interactions can happen. In our problem the dominant 
part of the LC propagator is DA = in a k^/[(k 2 +ie)(k-n+ie)], which gives rise to the noncausal 
interactions indicated by the l/(k ■ n + ie) factor from the second term of the propagator. So 
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when we speak of the choice of ie's we always refer to the second term of the propagator. Other 
components of the LC propagator are small. For example, in A-n = A + = gauge, D_ + (k) = 
and D (k) = 2ik_/k 2 k + is small compared to DA(k). 




(c) (d) 

Figure 1: Scalar particle production in light cone gauge involving two nucleons from each 
nucleus. The hard gluons come from nucleons 2 and 2', respectively. 

In order to illustrate the calculation we take two nucleons from each nucleus as shown 
in Fig. [TJ where the two nuclei are moving along the z-axis, and the dashed line represents 
the scalar field. As the two nuclei move towards each other at a very high energy, they are 
highly contracted in the direction of their motions. However, the nucleons inside a nucleus 
are well-ordered in the x_-direction, for the right-moving nucleus, or x + -direction, for the left- 
moving nucleus. We choose the nucleus 1 to be a right-mover which has a big plus momentum 
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component and is in A + = gauge. We label the front nucleon in the nucleus 1 to be nucleon 1 
which has a smaller ^..-coordinate than that of nucleon 2, xi_ < x^-. Similarly, nucleus 2 is a 
left-mover and is in A_ = gauge. The x + -coordinate of nucleon 1' is less than that of nucleon 
2', x' 1+ < x' 2+ . We will use (k + + ie) and (k- + ie) choices for the propagators in A + = and 
= gauge respectively, which means gluon fields propagate from the front to the back in a 
given nucleus. In other words, gluons in A + = 0(A_ = 0) gauge only propagate from a small 
x_ (x + )-coordinate to a large x_ (x + )-coordinate. With this very property a number of diagrams 
are prohibited at the beginning. A typical group of diagrams is shown in Fig. [TJ where we only 
draw one gluon from each nucleon. It is straightforward to add more gluons, but no more than 
two gluons from each nucleon in the quasi-classical approximation. If one wants to connect 



gluon k\ with, for example in Fig. 1(a) , gluon k, which come from nucleon 2 that has a smaller 
^-coordinate than that of nucleon 1', then one obtains a diagram that violates the (k + + ie) 
prescription. Such diagrams correspond to initial state interactions before the collision, they 
are suppressed by the mass M in the current ie prescriptions. A more detail discussion will 
be given in the appendix. Moreover, in the appendix we repeat the calculation in a totally 
different ie prescription, which leads to a huge amount of initial state interactions, then such 
diagrams are necessary in the calculation. In principle the current ie prescriptions would lead 
to severe final state interactions, however, color neutral scalar particle knows nothing about 
the gluon fields after the collision once produced, clearly it is not the case for gluon production. 
Therefore the four diagrams shown in Fig. [JJ are the only diagrams we have to calculate if we 
allow only one gluon from each nucleon. Moreover, hard gluons have to come from the last 
interacting nucleon in each nucleus according to the same reasons. In our approximation, we 
only need one hard gluon from each nucleus to make up the large mass M, and one-hard-gluon 
approximation is also one of the characteristics of the WW gluon fields. 

We will use the STW identities to evaluate the graphs in Fig. [TJ First, we focus on the left 



part of Fig. 1(a) The relevant factors are 



Pi + «e Pi+ + ie{k- Y + %e p 1+ - k + - ie , , 

u u 1 ^ ' 



9p,u 



ie K+ — ie 



k 2 + ie 

where V pj3ll = gf a bc[g/3 P (^Pi ~ *0a» + 0w*( 2 * ~ Pi)/9 + 9^{~k - Pi) J is the usual three-gluon 





(a) (b) 

Figure 2: Scalar particle production in the language of gauge rotation involving two nucleons 
(a) and an arbitrary number of nucleons (b) in each nucleus. A dashed line with an arrow 
denotes the gauge rotation. 

vertex, v v \ is the scalar-gluon vertex and rj-k = k + . The dominant part of the pi-propagator is 
proportional to + ie), while the (k — pi)-propagator contains (pi — &)p/G°i+ — k + — ie). 

The pi + -poles only come from the two propagators and are on opposite sides of the real axis. 
We now carry out an integration over p 1+ , the phase factor g^n-l^i--^-) w jth x\- — £2- < 
tells us to distort the contour in the lower-half plane and pick up the pole in the pi-propagator 
which sets p\ + = 0. Since pi_ is very small compared to I- for scalar particle production in a 
central rapidity region, we can replace p\^ by pip which allows us to apply the STW identities 
to the /5-index. The p\p factor will separate the vertex into two different parts 



PipF pPll = gfabc{ (g P ^k 2 - kpkp) - [g Pfl {pi - kf - {p x - fc) M (pi - k) p ] j. 



(4) 



The first term in Eq. (jlj) multiplied by the /c-propagator gives g pv —rj p k u / (k + +ie). However, 
the rj p term is eliminated by the {k — p^-propagator; i.e. {k — Pij^rjp = 0, so only g pv is left. 
This corresponds to eliminating the /c-propagator and bringing the pi-line to the gluon-scalar 
vertex. Similarly, the second term in Eq. 03]) eliminates the (k — p^-propagator and brings 
the pi-line back to the nucleon 2. Due to color factors the second term cancels with the graph 



shown in Fig. 1(b) Therefore, the non- vanishing contribution comes from the first term which 
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is denoted by a dashed line with an arrow, a notation introduced by t' Hooft [18] . attaching to 
the scalar-gluon vertex. Similar calculations can also be found in Ref . [TTHT3"] . The new pi-line 
changes the momentum of the (k — pi)-line to k, so there is no pi-dependence at the vertex 
v v \. This gluon line can be interpreted as a gauge rotation which is essentially a soft gluon, or 
a small- x gluon which is only coherent in the longitudinal direction. The gauge rotation brings 
color as well as transverse momentum but no longitudinal momentum to the vertex. We can 
also apply the same technique to the ki-\me which also becomes a gauge rotation. The sum of 
the four diagrams in Fig. [I] can be greatly simplified by the STW identities and become the 



one diagram shown in Fig. 2(a) with all the soft gluons interpreted as gauge rotations. Now 



the vertex part in Fig. 2(a) reads 



(* - Pi)„(l -k- h)iv vX (X (k — Pl )i{l -k- h)i [g uX k -(l-k)- k x (l - k) v ] 



(5) 



n-ik-pj-d-k-kjl+l-, 

where the transverse momentum components come from the two hard gluon propagators. We 
always assume that the mass squared of the scalar particle M 2 = I 2 is larger than any transverse 
momentum in the system. We see that although there are two soft gluons attached to the hard 
gluons, only the transverse momenta of the hard lines enter the vertex, which is the very 
characteristic of the WW field A^ w in momentum space. Therefore, we can identify the hard 
gluons from each nucleus as the WW fields at this order. This identification become even 
sharper by noticing that the gluon field emerging from one nucleus, either the left or right part 



of Fig. 2(a) , is the same as the result of Fig. 4 in Ref. 




Figure 3: Scalar particle production in nucleus- nucleus collisions. The gluon fields from each 
nucleus are interpreted as the WW fields. 
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This low order calculation can be easily generalized to an arbitrary number of interacting 
nucleons in each nucleus. The hard gluons still come from the last interacting nucleons while 
all the soft gluons are propagating backwards. We always apply the STW identities to the soft 
gluons that come from the most front nucleons first. The soft gluons become gauge rotations 



and merge at the scalar-gluon vertex as shown in Fig. 2(b) Adding more nucleons corresponds 
to attaching more gauge rotations to the hard gluons. All the gauge rotations can be resummed 
in coordinate space, then the hard gluons become the WW fields A^ w (x) of the nuclei [TT] . 
Therefore, one can first calculate the scattering amplitude at the lowest order where there is 
only one gluon from each nucleus as shown in Fig. [3J This gives 



M = -/(4p 1+ p 2 _)r i ^|f| 3 f ■ (6) 



We identify the gluon fields in Eq. ® as A^(k) = gT a k^/(k 2 k ± ) in A± = gauge p3], 
respectively. The cross section is given by 
da 9 2 2 2 f d 2 k 



9 i k + l - 1 ^MKW<(- k )M A i( l - k ) A i(~ l + *)]■ (7) 



d?l_dy 

Finally, going to coordinate space 

K(k,k + ) = J d 2 x,j d Xl . e i ^ +i ^- i ^A^(x 1 ,x 1 ^), 
Ai(k,k_)= [ d 2 x 2 [ dx 2 .e^ k - + ^-^A^(x 2 ,x 2+ ), 



one can replace A^ by the full WW AY[ w field which resums all the soft gluons. This gives 

= ^~tj d 2 xe u -^N la p(x)N 2af j(x) (a, = 1, 2), (9) 
which is a factorized form involving two unintegrated gluon distributions 

N ia ,(x) = ~J d\Tr[AVa W (h) A T(k + ^)] (< = 1.2), (10) 
where the index % labels the different nucleus. 

We will evaluate Eq. ( ITU]) following the procedure outlined in [T2] . In A + = gauge 
A^ w (x) can be written as 

A^ w (x,x_)= [ S{x,b„)T a S-\x,b_) ^^" p a {b, b_)9(b_-x^)d 2 bdb_, (11) 
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with 



S{x,x_) = Pexp!^igT a J In [\x - 6>]p°(6, 6_)0(6_ - x_)d 2 6d6_J, (12) 
as taken from [12]. Put A^ w in Eq. ffTUl) . and use 



(p a (6,6_)p a '(6:,^)) 



^=^(6 - 6>(6_ - fe'_)5 aa ,g 2 ^xG(x, Q 



(13) 



to average the color charges. In the McLerran-Venugopalan model we do not have either Q 2 - 
dependence or ^-dependence in Q 2 (d/dQ 2 )xG(x,Q 2 ). One finds 

S [ d % [d 2 b>db> ^-^ + ^ P ^^xGix^) 

(14) 



N a p{x) 



2tt 3 J -J ~ - \b-V\ 2 |6 + x-6'| 2 N* - 1 ^ ag 2 
x (Tr [5(6, 6 / _)T a S'- 1 (6, 6'_)5(6 + x, 6 / _)T a S'- 1 (6 + x, 6'_)]}. 



p(6', 6'_) = p re i is the normal nuclear density in the boosted frame and in light cone variables. 
One can do the d 2 lJ integration as 



d 2 b 



,(b - 6% (k + x — b[) 



\b — b' 1 2 1 6 + x — 6' I 



d 2 k e~ ik ^ 



k n k 



(k + v 



2\2 



-TV 



S a 0K o {\x\/i) + fi-Y-r-KtQxlfi) 



-7r[5 Q( g ln(|x|/x) + x a x ] 



(15) 



p is some infrared cut-off and x^ = x^/lxj is a unit vector. The trace is evaluated by expanding 
nucleons from 5(x, x_) one by one, the result is also taken from [12] 



(16) 



(Tr [5(6, b'_)T a S- l {b, b'_)S(b + x, b'_)T a S~\b + x, 6'_)] ) 

2 np Te iN c x 2 



CViV c exp 



-xG(x,l/x 2 )(6'_ +6' _) 



" 4(iV2-l)- 

where ±6q_ is the upper (lower) limit of the 6_ integration in Eq. f lT2|) . Plugging this back in 
Eq. (TH|) and performing the db_ integration and using the nuclear density in the center of mass 
frame p = p re i/7v / 2 with 7 the Lorentz contraction factor, one obtains 

iV 2 — 1 1 



Ni a p(x) 



An^aNrX 2 2 



<w + 



XaXp 

ln(|x|/i) 



d% (1 - e- Q te '*) 



(17) 



where Q 2 S are the corresponding saturation momenta for the two nuclei. For a spherical nucleus 



of radius R } it can be written as Q 2 = 8ir 2 aN c ^R 2 - b 2 pxG(x, l/x 2 )/(N 2 - 1). 
We can go to momentum space by 

NiapQu) = [ d 2 xe l ^ x -N ia p(x). (18) 
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The diagonal part of the above expression are the usual WW gluon distributions in momentum 
space, i.e. iVj(&J = SapN^^k^. The off-diagonal part are different gluon distributions, the 
so-called linearly polarized gluon distribution, also found in [20J 

nm = (2U* - Miwi.) = g£ J j^rjfk (i - ^")- (») 

fc ia are unit vectors in momentum space and J n (x) is the Bessel function of the first kind. 
We see that N a/3 (k) can be separated into two different unintegrated distribution functions, 
i.e. Niapikj) = ^SapNiikj) + (ki a ki/3 ~ | <W ) n i (i^ ) • This decomposition is usually encountered 
in studies of transverse momentum dependent factorization, see for example [2HI2H]- We can 
rewrite the cross section in terms of these two gluon distributions in momentum space as 

d 2 l_dy ~ 64 

A simple version of Eq. ( 120)) is found earlier by authors [19], where they calculate Higgs pro- 
duction in pA collisions. 

In the leading logarithmic approximation and when Z 2 ~ Q 2 s1 gluons with k 2 < Q 2 give the 
dominant contribution to the cross section and the factor ln(|x|/i) is large. We only need to 
keep the 5 a p term in Eq. (ITT)) , then obtain our main result 



J d 2 k 1 d 2 k 2 5 {2 \l_-k 1 -k 2 )^N 1 (k 1 )N 2 (k 2 ) + [2(fci ■ k 2 ) - ljm^iKfe 



(20) 



d 3 a gj 



2 7T 3 



d 2 ldy 16 

where 



d*xe**Nx(x)N2(x), (21) 



= J d% Tr[AV W (bJ ■ AY w {h + x)] = f d\ [I - e^W] , (22) 

are the unintegrated WW gluon distributions of the nuclei. It is not a surprise that WW gluon 
distribution manifest itself easily in this process, after all the WW gluon distribution counts 
the number of gluons in the nucleus and is the intrinsic gluon distribution of the nucleus. 

When I 2 > Q 2 , since the two nuclei contribute equally to the transverse momentum of 
the scalar particle, we also have fc 2 3> Q 2 . In this limit, |x| is very small, the exponential in 
Eq. (122)) can be expanded and the WW gluon distribution becomes additive 

Nl (k) ^ Ai^f^ = AN (k) (fc 2 »Q 2 ), (23) 
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where 

A = 2 J d 2 b lPi ^R*-b 2 (24) 
are the atomic numbers of the corresponding nucleus and 

N (k) = -42[xG(x,k 2 )], (25) 
7r dk 

is the unintegrated gluon distribution of a nucleon at the lowest order. The other gluon distri- 
bution has the same behavior, 



ni {k)^A^f^ (k 2 ^Q 2 s ), (26) 

7T rC 



where we have used / °° d\x\ \x\J2(\k\\x\) = 2/fc 2 . Therefore, in this limit we identify the two 
gluon distributions, then Eq. (120]) becomes 



d 3 



a 



2 

^ ! d 2 k 1 d 2 k 2 5 {2 \l_-k 1 -k 2 )2{k x ■k 2 )Ni{k 1 )N 2 {k 2 ), (27) 



d?l_dy 64 

where now A^(fc) take the form indicated in Eq. (|23|) . Similar expressions have also been ob- 
tained in a fcf-factorization formalism for pp [221123] and in the transverse-momentum-dependent 
factorization approach for pp [21] and pA [19]. 



4 Covariant gauge calculation 

Since we have used an unusual gauge choice to calculate the process, it is interesting to see 
what the calculation looks like in covariant gauge. As an example take two nucleons from each 
nucleus and use the same labeling as before. In covariant gauge interactions happen causally 
and locally. As the two nuclei approach each other, nucleon 1 meets nucleon 1' first. One way 
of producing the scalar particle is to have gluons from them creating the scalar particle, while 



nucleons 2 and 2' act as spectators as shown in Fig. 4(a) The hard gluons can also come from 



nucleons 1 and 2' as shown in Fig. 4(b) In this situation the gluon from nucleon 1 can be 
scattered by nucleon 1' before it meets the gluon from nucleon 2'. A similar process where the 
hard gluons come from nucleons 2 and 1' is also possible. It is not shown in Fig. |H Moreover, 



hard gluons can also be emitted from nucleons 2 and 2' as shown in Fig. 4(c) As the two 



nuclei pass through each other the gluon from nucleon 2 has to pass nucleon V first, and one 



13 



2 

o 



k-pi 



I - k 




(a) 



(b) 




k — p 



I — k — ki 




(c) (d) 
Figure 4: Scalar particle production in covariant gauge. 



(or two) gluon exchange can take place between them, that is the fci-line in Fig. 4(c) A similar 
process can happen to the other hard gluon, i.e. the pi-line. Finally, those two hard gluons 
collide and produce the scalar particle. This scattering sequence is equivalent to the diagram 



in Fig. 2(a) in light cone gauge. So, in covariant gauge, a hard gluon from a nucleon will be 



multiply scattered by nucleons from the other nucleus and receives large transverse momentum 
from the multiple scattering. In order to realize this picture calculationally it is important to 
include the phase factors which indicate which nucleons are in front and which are in back [16J. 
The factor is e l v\+(^\--x2-) e tki-[x 1+ -x 2+ ) w \ 1 [ c \ 1 dictates how we should distort the contours to 
pick up the poles in the propagators. 



Let us examine in detail the propagators of Fig. 4(c) In order to generate the large mass 
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M of the scalar particle the large plus momentum, l + , comes from nucleon 2 and the large 
minus momentum,/., comes from nucleon 2'. Momenta k\ and p\ are considered to be soft, 
so we may assume k-, \k\ <C I- and p\ + , \p \ <C l + . Moreover, we also have (k — pi)- ~ 
and (/ — k — ki) + w because minus(plus) momentum component of a gluon emitted from 
a right(left) moving nucleon is very small. Since k\ + ~ 0, we have k + / + . With such a 
momentum choice the picture of multiple scattering is fulfilled. The relevant integral is 

e ip 1+ (x!--X2-) e iki-(x' l+ -x' 2+ ) 



dpx + dkx- 



[(l — k — ki+ pi) 2 + ie] [(k —pi + h) 2 + ie] 

i \ J -i i i i \ (28) 

dPl+ n 1 1 ; WTTTi 1 — ~ / ^1 



4Z+Z_ J Pi+ -{l-k-k x + P X ) 2 /2L + 7 - (fc — p x + ^i) 2 / 2 ^+ + ^' 

where we have included two phase factors with xi_ — < and x' 1+ — x' 2+ < 0. With the same 
approximation other gluon propagators do not have additional poles contributing to Eq. (128j) . 
for example, we may take k\ + ie ~ — fc 2 + in the fci-propagator. The phase factors tell us to 
do both contour integrations J dk\_ and j dp\ + in the lower half plane to pick up the poles in 
the corresponding propagators. After the contour integrations, the two gluon propagators are 



put on shell, indicated by the vertical lines in Fig. 4(c) , then the two successive scatterings are 
independent. The k\- and pi-lines become soft and only carry transverse momenta and colors 
to the hard gluons which are similar to the gauge rotations we saw earlier in light cone gauge. 



Moreover, the phase factors also guarantee that Fig. 4(c) is the only non- vanishing diagram 



involving two nucleons from each nucleus. For example, with the same approximation used in 



obtaining Eq. ([28]) . diagram Fig. 4(d) gives 



dpi+dki-r, =-p- =- 

[(p 1 + k 1 ) 2 + ie][(l-k 1 -p 1 ) 2 + ie] 

/ e ipi + (xi--X2-) e iki-(x' 1+ -x' 2+ ) 
dkl ~ [2p 1+ V - (k, + p x ) 2 + ie] [2Z + (J_ - V) -(l-h- P,) 2 + ie] = °" 



(29) 



Since the fci_-pole in the second propagator lies opposite to the direction of the contour dis- 
tortion indicated by the phase factor e lkl ~^- x>1 +~ x ' 2 +\ Such diagrams give zero as required by 
causality in covariant gauge. 

It is straightforward to generalize the calculation to an arbitrary number of nucleons as 
shown in Fig. Suppose that the hard gluons are coming from zth and jth nucleons of nucleus 1 
and 2 respectively, the gluon from the zth nucleon in nucleus 1 can be scattered by the nucleons 
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Figure 5: In nucleus-nucleus collisions the hard gluons are multiply scattered and their trans- 
verse momenta are broadened. Gluons from nucleons i and j are effectively put on-shell after 
each scattering. 

coming before the jth nucleon in nucleus 2. Similarly, the gluon which comes from the jth 
nucleon in nucleus 2 can be scattered by the nucleons coming before the ith nucleon in nucleus 
1. One or two gluons exchanges can happen as the gluon passes each nucleon, and after each 
scattering the produced gluon is put on shell by a contour distortion. We see that a hard gluons 
has to travel through a certain length in the nuclear matter and its transverse momentum is 
broadened by the multiple scattering. This process can be described by a classical diffusion 
equation in momentum space [24J and the transverse momentum distribution can be found by 
solving this equation. The solution of gluon diffusing in nuclear matter is found in [12]. We 
will use the results derived in those two papers. The unintegrated gluon distribution at the 
final stage of multiple scattering can be written as 



where p is the nuclear density and is assumed to be uniform throughout the nucleus, N (x) is 



distribution for the gluon to have transverse coordinate x at a longitudinal position z, Zq is 
the longitudinal position where the hard gluon finally emerges from the nuclear medium and is 
ready to produce the scalar particle. For the gluon from ith nucleon we can take zq « Zj-x- In 




(30) 



the initial gluon distribution which can be taken to be xG(x, 1/x 2 ), f(z,x) is the probability 
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the case of gluon p r broadening, for a spherical nucleus, f{z,x) can be written as 



z + VR 2 - b 2 2 , 



/(*,g) = exp|- 8 ^ 2 _ &2 g'ffj, (31) 

where is the radius of the nucleus and 6 is the impact parameter. Since the hard gluon can 
come from any nucleon in the nucleus, we have to sum up all the possibilities which corresponds 
to the dzo integration in Eq. (|30|) . It gives us 



ivi(xj oc / ctz exp < -x Q s 



VR 2 -b 2 



8VR 2 - b 2 



k2 ^? (1 -^ /4) ' 



(32) 



which is the WW type gluon distribution and is the same as Eq. (I22j) . This one-to-one corre- 
spondence between the LC calculation and the covariant calculation is also explained in [T2]. 
It is now quite convincing that WW gluon distribution is the right gluon distribution that 
should be used for this process. It is interesting to note that in covariant gauge, as in light 
cone gauge, the nucleons "behind" the nucleon which gives the hard gluon creating the scalar 
particle are viewed as non-interacting, which is the very property that gives rise to the WW 
gluon distribution. However, in covariant gauge, the diagrams do not clearly indicate that the 
cross section is factorizable. 
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(a) (b) (c) 

i i 
i i 




(d) (e) 

Figure 6: Collisions are viewed in a different choice of ze's. The hard gluons come from nucleons 
1 and 1'. The soft gluons from nucleons 2 and 2' become the gauge rotations. 

Appendices 



A Other choices of ze's 



In this appendix, we will show how a different regularization in the LC propagator gives the 
same result. Although there are three more cases in addition to the one discussed in the body 
of the paper, here we will only study the most complicated one where the gluons in the nuclei 
propagate in the negative x_(x + )-direction in A + = 0(A_ = 0) gauge. The two remaining cases 
can also be shown to give the same result by applying the same technique. Now the choices of 
ie's are taken to be (k + — ie) and (k- — ie) for A + = and A_ = gauge respectively, which 



18 



means that gluons in A + = 0(A_ = 0) gauge propagate from a large x_(x + )-coordinate to a 
small x_(x + )-coordinate. With all the gluons propagating in this way initial state interactions 
are complex and entanglement between the soft gluons can occur. Even nucleons belonging 
to different nuclei could affect each other, therefore the factorizability is not clear at the very 
beginning. This is quite different from what we have seen previously where there is no initial 
state interactions between the two nuclei. We will see that the STW identities guarantee that 
the cross section is still factorizable and is the same as the one obtained previously. 




(d) (e) 

Figure 7: Additional diagrams that are needed in the case where the hard gluons are from 
nucleons 1 and V. 

In light cone gauge we have the privilege to specify where the hard gluons come from, so 
in order to simplify the calculation we further classify the diagrams by how the hard gluons 
are connected to the nucleons and group them into several distinct classes. The advantage of 
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doing so is that, on one hand, one can avoid double counting diagrams and, on the other hand, 
the gluons that we have to add to the diagrams should be soft so that we can apply the STW 
identities to them. 

We start with the case where the hard gluons are coming from the first nucleons on each 
side as shown in Fig. |6j There are still two more soft gluon lines that should be added to the 
graph. In order to avoid double counting we can fix one of the soft lines first, for example, 



the gluon line connecting nucleons 1 and 2 in Fig. 6(a) and Fig. 6(b) Then draw different 



connections of the second soft line so that we can apply the STW identities. If we add Fig. 6 (a 



and Fig. 6(b) together the soft gluon from nucleon 2' becomes the gauge rotation shown in 



Fig. 6(c) Then we focus on the first soft line and enumerate its possible connections, now as 



a dashed gauge line, while doing the same thing to the second soft gluon. We obtain another 



diagram shown in Fig. 6(d) Adding Fig. 6(c) and Fig. 6(d) together, we change the soft gluon 



from nucleon 2 to a gauge rotation and obtain the diagram shown in Fig. 6(e) Since the 
gauge rotation acts independently on the hard gluons, the gauge rotation from nucleon 2' can 
also attach to the hard gluon from nucleon 1. One obtains an additional diagram shown in 



Fig. 7(a) We need one more diagram, Fig. 7(b) , in order to apply the STW identities. But 



such a diagram with a three-gluon-scalar vertex is suppressed by a power of mass squared of 
the scalar particle M 2 = I 2 and can be neglected in our approximation. It is shown here only 



for the sake of completeness. If one add all the diagrams Fig. 6(e) , Fig. 7(a) , Fig. 7(b) together, 
the STW identities guarantee that they give zero contribution. A similar argument can also be 



made to obtain the diagrams shown in Fig. 7(c) , Fig. 7(d) and Fig. 7(e) , their sum is also zero. 



Therefore the soft gluons from the last interacting nucleons do not contribute to the collisions 
in this specific choices of ze's. This result is what one should expect from the STW identities. 
Since the identities tell us that if we enumerate all possible insertions of a longitudinal polarized 
gluon line to a certain graph their sum should be zero. This is exactly what we have in this 
case, where the gluon fields coming from either nucleon 2 or 2', the last interacting nucleon, can 
gauge rotate everything that comes before them. Therefore all possible connections should add 
up to zero. So it is not surprising that even though we might have a huge number of diagrams, 
gauge invariance guarantees that many of them add to give no contribution. 
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(a) 



(b) 



(c) 






(d) (e) (f) 

Figure 8: Other two different cases where the hard gluons come from nucleons 2 and 1' and 
nucleons 2 and 2'. The later one gives the dominant contribution. 

Similarly, one can also consider the case where the hard gluons are coming from nucleons 



2 and 1', shown in Fig. 8(a), Fig. 8(b) and Fig. 8(c), and the symmetric case, where the hard 



lines are coming from nucleons 2' and 1, which is not shown in the figure. Since these two cases 
are essentially the same, it is sufficient to study only one of them. Due to the STW identities, 
they also give us no contribution. 

Finally, the last case is the one where the hard lines are coming from nucleons 2 and 2'. 
Since now the soft gluons can gauge rotate only certain parts of the diagram, it should give a 
nonzero result. If there is an entanglement between the two soft gluons, one obtains diagrams 



like Fig. 8(d) and Fig. 8(e) However, this kind of diagrams are suppressed by M 2 . Since the 
three-gluon vertex brings in an additional factor of transverse-momentum squared, a factor of 
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M 2 comes from the denominator to make up the right dimension. Therefore those diagrams 
have small contributions in our approximation. Moreover, one could also make a contour 
distortion explained in the appendix to make such diagrams zero. If the two soft gluons are 



directly connected to the hard gluons, one obtains Fig. 8(f) which is the dominant contribution 
in this choice of ie's. One can immediately recognize that this diagram is exactly the same as 



the one in Fig. 2(a) 



B Suppression of the initial state interactions 



k-pi 



I — k — k\ 




Figure 9: The initial state interactions are suppressed by the heavy mass of the scalar particle. 

In Section [3] we argued that the initial state interactions are suppressed by the scalar mass 
M in the (k± + + ie) and (k- + ie) gauge choice. Here we will do an explicit calculation to show 
that it is true. We take the example mentioned in Section |3] The diagram is shown in Fig. |9l 



which is quite similar to Fig. 1(a) but the connection of ki-\ine is changed. Let us focus on the 



part of the diagram that involves the gluons from nucleus 2, nucleons 1' and 2'. The relevant 
factors in the integrand are 



r] a k 
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k{ + ie k\- + ie 
i 



9v» 



T]fj,(k + k 



1 v 



r] v (k + k 



(k + h) 2 + ie |_ (k + k 1 ) + + ie (k + k l ) + -ie_ 
f}„(l -k- fci)£ 



(33) 



(/ - k - h) 2 + ie {l-k-k x )- + ie 
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where fj-k = k_. Since the big momentum component that flows through the (/c+/ci)-propagator 
is the big plus momentum component coming from nucleon 2, it is natural to use A + = gauge 
for the (k + fci)-propagator. We can distort the contour to the upper half /ci_-plane to pick up 
the k\- = (I — k)- + ie pole. Since pi_ ~ and (k — ~ 0, then (I — A;)_ w Z_ ~ M is large 
in the center of mass system of the colliding nuclei. In order to compensate for the + ie 
and (A; + ki) 2 + ze ~ 2k + k\- — (k + fc 1 ) 2 + denominators, we would need a coming from 
the numerator, a k\_ cancels the two denominators and one k\_ is necessarily to form the large 
mass M. There are only two vertices, and v u \, that involve However, the v v \ vertex 
cannot give a factor of k\- because of (I — k — ki)j^. Therefore, there is sufficient convergence 
in the fc^-plane for the contour to vanish at infinity, and such diagrams are suppressed by the 
mass of the scalar particle. We could arrive at an even stronger result more directly by noticing 
that we can replace (I — k — by (/ — k — kx)\ after the contour distortion and then use 
the current conservation in the A-index to obtain a zero result. However, such a diagram is not 
suppressed in — ie) choice for the propagators in A_ = gauge. This new ie prescription 
changes the second and the last factors in Eq. (!33l to — ie) and 1/[(Z — k — fci)_ — ie], 

respectively. We can still distort the contour in the upper half plane to pick up the k\- = ie 
pole. All the other /c!_-poles lie on the other side of the real axis, and no suppression occurs in 
this situation. After the contour integration, one can replace k^g by kip and then use the STW 
identities to simplify the diagrams. 
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